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Abstract 

The overcompleteness of the coherent states basis leads to a multiphcity of representations of 
Feynman's path integral. These different representations, although equivalent quantum mechani- 
cally, lead to different semiclassical limits. Two such semiclassical formulas were derived in for 
the two corresponding path integral forms suggested by Klauder and Skagerstan in 2]. Each of 
these formulas involve trajectories governed by a different classical representation of the Hamilto- 
nian operator: the P representation in one case and the Q representation in other. In this paper 
we construct a third representation of the path integral whose semiclassical limit involves directly 
the Weyl representation of the Hamiltonian operator, i.e., the classical Hamiltonian itself. 

PACS numbers: 03.65.Db, 03.65.Sq 
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I. INTRODUCTION 



In recent years there has been a renewed interest in semiclassical approximations with 
coherent states. These approximations have found apphcations in many areas of physics and 
chemistry. The semiclassical coherent state propagator has a long history, that starts with 



Klauder 0,14,0] and Weissman |6[. Several properties of the propagator were su 



studiec. 



for a number of fundamental quantum processes (see, e.g.. 
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1^). More recently, a detaUed der.vat.on of the sem.class.cal propagator for systems 
with one degree of freedom was presented in |1||. 

The set of coherent states forms an non-orthogonal over-complete basis, since each state 
in the set can be written as a linear combination of the others. This overcompleteness, on 
the other hand, has important consequences in the path integral formulation of the propaga- 
tor. It implies in the existence of several forms of the path integrals, all equivalent quantum 
mechanically, but each leading to a slightly different semiclassical limit. Klauder and Skager- 
stam proposed two basic forms for the coherent state path integral, each of them having 
their corresponding advantages and problems j^. The semiclassical limit of these two basic 
propagators were considered in 0] where it was shown that both propagators can written in 
terms of classical complex trajectories, each governed by different classical representation of 
the Hamiltonian operator H: the P representation in one case and the Q representation in 
other. We briefly review these representations in section 2. The phase appearing in these 
semiclassical formulas turns out to be not just the action of the corresponding complex clas- 
sical trajectory, but it also contains a 'correction term' I that comes with different signs in 
each formula (see Eqs. (fT3j) and (fTBj)). 

In it was also suggested that a semiclassical representation involving directly the Weyl 
representation of H, or the classical Hamiltonian, could probably be constructed, and a 
formula for this representation was conjectured. In this paper we derive this formula by 
constructing a new representation of the quantum mechanical path integral, Eq. ()2(j|l . and 
deriving its semiclassical limit Eq.()5H). We show that the classical trajectories involved 
in this formula are governed by the average between the P and Q representations of the 
Hamiltonian operator. The correction term in the phase, on the other hand, turns out to 
be one half of the difference between the corresponding terms in the previous formulations. 
We then show that this average Hamiltonian can be replaced by the classical one and the 
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correction term discarded, the error being of order h"^. Our final result is the conjectured 
formula, Eq.(jnZj). 

This paper is organized as follows: in section 2 we review the path integral constructions 
of Klauder and Skagerstan and their semiclassical approximations 3]- In section 3 we 
construct the new quantum representation and derive its semiclassical limit. Finally in 
section 4 we show that, within the validity of the approximations, this formula coincides 
with the equation suggested in [jj]. 

II. THE COHERENT STATE PROPAGATOR AND ITS SEMICLASSICAL AP- 
PROXIMATIONS 

A. The propagator 

The coherent state 1^;) of a harmonic oscillator of mass m and frequency uj is defined by 

|z) =e-5l-l'e"'''|0) (1) 
with |0) the harmonic oscillator ground state and 

In the above q, p, and are operators; q and p are real numbers; z is complex. The param- 
i_ 1 

eters b = {h/muj)^ and c = (hmu)^ define the length and momentum scales, respectively, 
and their product is h. 

For a time-independent Hamiltonian operator H, the propagator in the coherent states 
representation is the matrix element of the evolution operator between the states \z') and 
\z"): 

K{z",T;z',0) = (/'|e-i^^|z')- (3) 

We restrict ourselves to Hamiltonians that can be expanded in a power series of the creation 
and annihilator operators and d. 

In the derivation the semiclassical limit of the propagator, the Hamiltonian operator H is 
somehow replaced by a classical Hamiltonian function H{q,p). This 'replacement', however, 
is not uniquely defined, and the ambiguities that exist in the relation between the operator H 
and the function H{q,p) also arise in connection with the overcompleteness of the coherent 
state basis, as we shall see in the next subsections. 



There are actually many ways to associate a classical function of position and momentum 
A{q, p) to a quantum mechanical operator A. However, three of them are specially important. 
The first one, denoted AQ{q,p), is called the Q representation of the operator A and is 
constructed as follows: one writes A in terms of the creation and annihilation operators 
and d in such a way that all the creation operators appear to the left of the annihilation 
operators, making each monomial of A look like Cnmot^a™. Then we replace d hj z and 
at by z*. The inverse of this operation, that associates a quantum operator to a classical 
function, is called 'normal ordering'. In this case one first writes the classical function in 
terms of z and z*, with all the z* 's to the left of the z's, and then replace z by a and z* by 

at. 

The second possibility, called the P representation of A, is obtained by a similar procedure, 
but this time the monomials of A are written in the opposite order, such that they look like 
Cnm^f^Q^t™-. Once the operator has been put in this form one replaces again a by 2; and at by 
z* to obtain Ap{q,p). The inverse of this operation is called 'anti-normal ordering'. Notice 
that the differences between the two representations come from the commutator of q and p, 
which is proportional to h. Therefore, these differences go to zero as h goes to zero. 

There is, finally, a third representation which is the most symmetric of all, and therefore 
the most natural. It is given by the Wigner transformation 



s 



Aw{q,p)= / dsenP^{q-- A q + -) . (4) 
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Aw{q,p) is called the Weyl representation of A Its inverse transformation consists 

in writing the classical function in terms of z and z* considering all possible orderings for 
each monomial and making a symmetric average between all possibilities before replacing z 
and z* by the corresponding operators. As an illustration of these three representations we 
take 

1 1 2 4 

^=2^ + r +^ 

{m = h = 1) for which we obtain 

Hq = + +x^ + \{P + 6-2) + 3PX^ + 3574 

Hp = + x^) + - + r^) - 36^x2 + 3574 
where b is the width of the coherent state. Notice the term proportional to x"^ that appears 



with opposite signs in Hq and Hp, really modifying the classical dynamics with respect to 
Hw 

In the next subsections we shall see how these different representations appear naturally 
in the semiclassical limit of the coherent state propagator. 



B. Basic Path Integrals and their Semiclassical Approximations 

The calculation of the semiclassical propagator in the coherent state representation start- 
ing from path integrals was discussed in detail in In this section we summarize these 
previous results emphasizing the non-uniqueness of the semiclassical limit as a consequence 
of the overcompleteness of the coherent state representation. The reader is referred to 
for the details. 

In order to write a path integral for K{z", T; z', 0), the time interval has to be divided into 
a large number of slices and, for each slice, an infinitesimal propagator has to be calculated. 
As pointed out by Klauder and Skagerstam Q], there are at least two different ways to do 
that. Each of these gives rise to a different representation of the path integral. Although they 
correspond to identical quantum mechanical quantities, their semiclassical approximations 
are different. We review the construction of these two representations below. 

The first form of path integral is constructed by breaking the time interval T into 
parts of size r and inserting the unit operator 

f , , d^z , , 

t= / \z) — {z\ (5) 



TT 

everywhere between adjacent propagation steps. We denote the real and imaginary parts of 
z hj X and y, respectively. In all integrations, d'^z/ir means dxdy/71. After the insertions, 
the propagator becomes a 2(A^ — l)-fold integral over the whole phase space 

„ N~l ,2 ^-1 

Kiz",t;z',o) = / {n ^} n{(-.>ii^"^^^*^^i-.)} (6) 

with zn = z" and zo = z'. Using the coherent state overlap formula 

{zj+i\zj) = exp + z*^^Zj - ^kjf | (7) 

and expanding e~^^'^^^ ^ 1 — iHr /h we write 

f 1 1 1 

{Zj+i |e-^^(*^)"|2;,-) = exp 1 2 " 4)^^' " 24+i(^J-+i - z,) - y^i+u j (8) 



where 

{zj ^i\H{t,j)\z j) 

{Zj+i\Zj) 

and {l — iHj+ijT/fi) has been approximated again by e~*^J+i'-'^/'*. With these manipulations 
the first form of the propagator, that we shall call Ki, becomes 



Mz",t;z',0)= /{n^}-P E 



N-l ,a fN-1 

1 



d^ 

j=l k i=0 



(10) 



When the limit N ^ oo (respectively r — 0) is taken, the above summations turn into 
integrals, and expressions (jlOj) appears to be exact, were it not for the well-known problems 
attached to the meaning of such functional integrals. Also, turns into the smooth 

Hamiltonian function T-C{z,z*) = Hi{z,z*) = {z\H\z). Using the properties a\z) = z\z) and 
{z\a) = {z\z* , we see that TC can be easily calculated if H is written in terms of creation 
and annihilation operators with all 's to the left of the d 's. Therefore, Ti is exactly the Q 
symbol of the Hamiltonian operator jlS^ . 

The second form of path integral starts from the "diagonal representation" of the hamil- 
tonian operator, namely 

H = [ \z)H2{z\z) — {z\. (11) 



TT 

Assuming that H is either a polynomial in p and g or a converging sequence of such polyno- 
mials, this diagonal representation always exists. The calculation of H2 is not as direct as 
that of Hi, but it can be shown 18] that H2{z*, z) is exactly the P symbol of H. This second 
form will be contrasted with the first-form hamiltonian function Hi{z*, z). To facilitate the 
comparison between the second form of the propagator, that we call K2, and the first form 
Ki, it is convenient to break the time interval T into — 1 intervals, rather than A^. We 
write 

N-l 

K2{z",T-z',0) = {z"\l[e-f^\z') (12) 

i=i 

and, following Klauder and Skagerstam, we write the infinitesimal propagators as 

e-^-^ - / \z,) (1 - jH2{zl z,)) - / k,)e-x-^(^^^^)^(z,| . (13) 
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The complete propagator K2 becomes 

f^~^ d^z- {it 1 



i=i 

j=l I j=0 



1 1 ir ' 

2 (4+1 ~ ~ 2^j+i(^j+'^ " ^3) " -^H2{Zj, Zj) 



(14) 

The differences between Ki and K2 are subtle but important. While the two arguments 
of Hi in Ki belong to two adjacent times in the mesh, the two arguments of H2 in K2 belong 
to the same time. Although both forms should give identical results when computed exactly, 
the differences between the two are important for the stationary exponent approximation, 
resulting in different semiclassical propagators. The semiclassical evaluation of Ki and K2 
were presented in detail in 111. Here we only list the results: 

Mz",t;z',0) = exp + - i(|/f + |.f )| , (15) 

KAz", t, z',0) = J: ^i^, exp U{S,. - I,.) - l(\z'f + \zf) ] , (16) 

where 

S^^ = S^v", u', t) =J dt' [^{iiv - vu) - H,{u, V, f)] - f{u"v" + u'v') (17) 



is the action and 

-dt (18) 



2 Jq dudv 

is a correction to the action. The sum over u represents the sum over all (complex) classical 
trajectories satisfying Hamilton's equations 

inu = 



dv 
OH 

ihv = — (19) 

ou 

with boundary conditions 

u{0) = z' = u' , v{t) = z'"- = v" . (20) 

The factors Jj are an important part of the above formulas and they are absolutely necessary 
to recover the exact propagator for quadratic Hamiltonians. If one neglects it, even the 



Harmonic osci 
see refs. 



m 



lator comes out wrong. For a discussion about non-contributing trajectories. 



C. The Conjectured Weyl approximation 



As discussed at the begining of this section, for a given quantum operator the first- 
form Hamiltonian is given by Hi{z*, z) = {z\H\z). It can be obtained by writing H in terms 
of the operators d and in normal order, so that each monomial in H look like CnmS^"a™. 
Replacing by and d by z yields Hi. H2 can be calculated by writing H in anti-normal 
order, where now each monomial looks like c„mS"d^™', and then replacing d^ by z* and d 
by z. A third type of Hamiltonian function can be obtained from H by using the Wigner 
t r ansf ormation : 



s 



s 



Hw{q,p) = J dsen^^i^q-- H q + -) . (21) 

This is the Weyl Hamiltonian. Since H]^ is obtained from H by completely symmetrizing 
the creation and annihilation operators, it turns out that H]^ is an exact average between 
Hi and H2 if H contains up to cubic monomials in d and d^, but only an approximate 
average for other cases. The semiclassical formula with Hi comes with a correction +Ji to 
the action and that with H2 comes with a correction of —l2- This suggests a third type of 
semiclassical approximation for the propagator, where one uses the Weyl Hamiltonian and 
no correction term, since the average of +/i and —Jo should be approximately zero. This is 
the Weyl approximation, which was conjecture in jlf: 



with Sw given by Eq. lfTTj) with Hi replaced by Hw- 

Of the three semiclassical approximations presented, the Weyl approximation seems to 
be the most natural, since it involves the classical hamiltonian directly and no corrections 
to the action. However, this formula does not follow from the two most natural forms of 
path integral proposed by Klauder and used in this section. In the next section we propose 
a third form of path integral whose semiclassical limit is indeed the Weyl approximation. 

III. A MIXED FORM FOR THE PATH INTEGRAL 

The new form of path integral we describe in this section is based on the fact that Hw is 
almost the average of Hi and H2. The idea is to force this average to appear by combining 



8 



the first and second form of path integrals in alternating time steps. We start from 

N-l 

K{z",t;z',0) = {zr,\ J] e-^^^^>o) (23) 

j=0 

where = z" , zq = z', Tj is the time step and we take to be even for convenience. 
Although we shall consider the time steps tj to be all equal later, we keep the index j for 
now to keep track of the time intervals. 
For j odd we approximate 



e-i"^^ - j [l - fH,iz],z,)j ^(.,1 ^ j |.,)e-^^^(^l-^)^(.,| . (24) 
For j > even we simply insert a unit operator on the right of the infinitesimal propagator: 

e-i^^=|e"i^^|^)^(.,|. (25) 

Multiplying this operator on the left by the bra {zj-i\ coming from the odd term j — 1 and 
using the approximation employed in the first form of path integrals, Eq.(jHl), we get the 
following mixed form for the propagator: 

K{z^,T-zo,o)= y exp x: 



IT IT 

■—ajH2j - -^bjHij 

N-l 



where aj is zero for j even and one for j odd, bj is zero for j odd and one for j is even. At 
this point we have suppressed the index on the time intervals and have taken tj = t. The 
exponent / is given by 

N-l . ^ ^ ■ 

f{z\ z) = J2\ 2 (4+1 - 4)% - 24+i(^.+i - z,) - ja,H2, - jb,H,,, (27) 
where we have introduced the abbreviated notation = H2{zj, zj) and Hij = Ti.{z*_^^^, Zj). 



A. The Stationary Exponent Approximation 

In the semiclassical limit ^ — > we can approximate the integrals ()26|) by looking for the 
places where the exponent / is stationary and replacing it in their vicinity by a quadratic 



form of its variables {z*, z). We find the stationary points by requiring the vanishing of the 
derivatives of / with respect to z and z* separately. We obtain 

dz,-""^^^ n dz, h dzj •^-^'•••'^^ ^ 

df iajj^iT dH2j ibjrdHi j 



We now introduce new integration variables t] and rj*, which describe the deviations from 
the points of stationary exponent, z —>■ z + t] , z* ^ z* + ri*, with the boundary conditions 

Vo = Vo = m = VN = . (29) 

Expanding the exponent into a Taylor series in {ri*,ri) around the stationary points {z*,z) 
up to second order and re-inserting the result into (jHl) yields 



The integrals in Eq. ()30|) can be carried out using the same techniques presented in Q|. 
The idea is to integrate first over 7]\ and r/i, then over r/2 and 772? etc. A recursion formula 
can be readily established and, once all integrations are done, we obtain 

N-l 

KizN, t; zo, 0) = e^(^*'^) J] , = (31) 



where Xj satisfies 



' 2% dzf 2h dzf ^ ' 



V h dz*dzj_i J 



^ h dzjdzj_i^ 

h dz^^,dz,J ^V""'-'^^^^^-'^^^^'-' 

for J = 1, . . . , AT - 1 with Xq = 0. 
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B. The Effective Hamiltonian 



To obtain the continuum limit of the discrete equations of motion (|28|) we first note that, 
since aj is zero for j even and one for j odd, and hj is zero for j odd and one for j is even, 
the second of Eqs. gives 

-Zj+i + Zj = {) (34) 

for j even. 

This result motivates the choice of a new time step e = 2r, in such a way that the 
discretized time evolution goes from Zj directly to Zjj^2 for j even and from 2;* directly to z*^2 
for j odd . This choice of time step makes sense if one can find an effective Hamiltonian able 
to perform the corresponding evolution for both variables z and z* . Indeed, the equations 
of motion can be put the form 



Zjj^i — Zj-i i dH, 



e h dZj 

4+1 ~ 4-1 _ ^ '^Hefj 



h dzj 



l,3,5...,A^-3 

2,4,...,iV-2. (35) 



where 



^efd = (^i-^ij + bj-iHij_i + aj^iH2j + ajH2j+i) /2 

Hi, + H,, . (36) 

tor ? even 

2 

H1J-1 + H2J+1 . -I-, 

— — for ? odd . 

2 

Note that, because we are skipping points with the new time step e = 2r, we miss the point 
z^ in Eq. ()35|) . This, however, does not affect the limit of the continuum, as long as we take 
= z'^ to ensure the proper boundary condition. 
In the limit where 2r = e goes to zero the effective Hamiltonian reduces to 

n,(.',.)^^h(£A±I^2l£i£l, (37) 

As in the case of the semiclassical formulas Eqs. ()15|) and (fT^ for Ki and K2, the station- 
ary trajectory is usually complex. It is therefore convenient to follow the notation introduced 
in [l| and make the substitutions 

1 /q .p\ 
^/2\b c 



(38) 
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In terms of u and v, the stationary phase conditions ()28|) turn into Hamiltons equations 
with He- 

dHc 



ihit = + 



dv 



ihij = -^ (39) 
ou 

with boundary conditions identical to Eq.((2nj. The function / can also be simplified to 



f= dt' 



* 1 i 

-iyu - iiv) - -^^c{u, V, t') 



+ -{v"u" + v'u')--{\z"Y+\z'f) (40) 



where m(0) = u' = z\ f (0) = v' , u{T) = u" and v{T) = v" = z"*. 

Next we calculate the product appearing in (jHlll . Performing the limit N ^ oo and using 
the expansion ln(l + x) = x + O(x^) we obtain 



j=l 3 1 3 j 

In order to transform these sums into integrals we note that 

hm Yl ^i^^Fih)) = lim r(F(tO + Fit,) + Fit,) + ...) 

N—>-oa ' N^oo 

^ j=l ^ 

= \ lim r(F(ti) + F(ti) + ^(ta) + Fit,) + ^(tg) + -..) 

= ]- lim r(F(ti) + F(t2) + Ffe) + Fit^) + ^(tg) + •.•) 

t 

= \jdt'Fit') (42) 



since, for smooth functions, -F(tj) — * F(tj+i) as r — > 0. The integrals with the coefficients 
hj also acquire the 1/2 factor. Using these results we obtain 

r = = expj-- J di' (-—(*') + -(— + ^)A-(t'))| 
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C. The Effective Phase 



Replacing T and / into Eq.(jnH) we obtain 



K{z",t-z',0)=exp(^~ I dt' 



it') + ^(t')xit')] 



dudv 



^ {vu - uv) - ^Hc 



+ liv'u' + v"u")-^{\zf + \zr)\. (44) 



exp| / dt' 



We still have to write the continuous form of the discrete recursion formula ()32|1 for X{t) 
In the limit — oo we obtain the nonlinear differential equation 



Xit) 



i d^Hc 



(45) 



2h dv"^ h dudv 
with the initial condition X{0) = 0. 

This equation was solved in [l| and the result is X = where 5u and 5v are solutions 
of the linearized Hamilton's equations 



6u 
6v 



h dudv 
i d^Hc 



5u — — 



5u+ - 



h dv"^ 
i d^Hc 



6v 



6v 



(46) 



fi du^ ' h dudv 

where the derivatives are calculated at the stationary trajectory and the initial conditions 
are Su{0) = and 6v{0) arbitrary. The second term in the first exponential of ()44|) can be 
now transformed with the help of 



i d^Hc 



X 



i d'^Hc 6u 16v i d^Hc 1 d 



h du^ 2h du^ Sv 2 6v 2h dudv 
so that the first exponent of ()44|1 becomes 

I f'^^,^ld''H2,^„ , d'Hc,^,, 
'o 

i 



-—\nSv 
2dt 



i d^Hc 
2h dudv 



exp 



exp < + 



2 dudv 



du 



(47) 



(48) 



d^Ho 



dudv 



it')] exp 



dt' 



dt' 



-{ln6v) 



Sv' 



6v' 



-exp 



h 



Sv' 

17' 



exp 



h 



(49) 



The pre-factor Sv' jbv" can also be written in terms of the action using dSc/du' = —ihv', 
where 

t 



Sc{v",u',t): =Jdt'[f{u 





uv — vu] 



Hciu,v,t' 



'}{u"v" 



u'v') 



(50) 



13 



is the effective action. In the end we obtain 
where 

Ic = \{h - h) (52) 

is the effective phase. 



IV. THE WEYL APPROXIMATION 

Equation (jKT|) is the main result of this paper. It represents a third semiclassical 
approximation for the coherent states propagator, involving the effective Hamiltonian 
He = {Hi + H2)/2 and the effective phase Ic = (h - h)/'^. 

For the harmonic oscillator Ji and I2 are exactly equal and Iq = 0- In this case He coin- 
cides with the classical, or Weyl, Hamiltonian Hw, and the conjectured Weyl approximation 
(I22|l is obtained. In fact Ji = I2 and He = Hw for all polynomial Harniltonians involving 
up to cubic powers of g or p [l^ . This can be seen from the formulas 



Hiiz\z)=exp[h]Hw{z*,z) (53) 



H2{z\ z) = exp ( -1-6 ) Hw{z\ z) (54) 



2 



where 5 = d /dz*dz. This gives 



He = cosh (^Is^Hw = Hw- Ig^^Hw + ... (55) 

which shows explicitly that He = Hw for up to cubic polynomials. Besides, using the 
relations 

q={z + zn{^) P=(^-^*)(^) (56) 

with 6 ~ c ~ 0{h^^'^), we see that quartic or higher order terms contribute to He or le 
only terms of order h^. These terms can in principle be neglected, since they are beyond the 
scope of the approximation. With these considerations we can rewrite the propagator (plj) 

as 



i d^Sw fir, 1/1 /|2 , I //|2^ 
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where Sw is given by Eq. (fSUj) with He replaced by Hw. This is the Weyl formula conjectured 
in 

As a final remark we notice that the differences between the three semiclassical formulas 
presented in this paper are of the order of and go to zero in the semiclassical limit. 
These differences, however, are always relevant at low energies, and can be made explicit 
by considering the Fourier transform of these time dependent formulas. A discussion of 
the energy representation of these semiclassical propagators were presented in section 6 of 
ref.jll, including the derivation of semiclassical quantization rules, and we refer to it for 
further details. 
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